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Abstract 

— |- We implement a direct detection scheme based on hybrid photodetectors to experimentally in- 

o 

vestigate high-order correlations for detected photons by means of quantities that can be experi- 
mentally accessed. We show their usefulness in fully characterizing a multimode twin-beam state 
in comparison with classical states and, in particular, we introduce a nonclassicality criterion based 

> 

on a simple linear combination of high-order correlation functions. Our scheme is self-consistent, 
allowing the estimation of all the involved parameters (quantum efficiency, number of modes and 
average energy) directly from the same experimental data. Results are in very good agreement 
with theory, thus suggesting the exploitation of our scheme for reliable state characterization in 
quantum technology. 
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Correlations play a fundamental role in the investigation of optical coherence pQ. Since 
their introduction [2] , correlation functions have been extensively investigated in connection 
to quantum state characterization, to define nonclassicality criteria [3J, and for the enhance- 
ment of ghost-imaging protocols pEJ [5] . From the experimental point of view, the study of 
this topic dates back to the pioneering work in which Hanbury Brown and Twiss discov- 
ered photon bunching in light emitted by a chaotic source [6]. In the last decade, many 
experiments have been developed in which the photon-number correlations have been used 
to characterize the entangled states generated by parametric down- conversion [7HTU]. In all 
the mentioned cases, the detection was performed by means of single or arrays of avalanche 
photodiodes so that the possibility to recover the correlation of the number of photons was 
quite straightforward [HI [12], although the intensity range actually investigated by these 
systems was limited to much less than one mean photon [13]. However, the investigation of 
more intense light beams is of extreme interest, especially in view of possible applications to 
quantum technology: pulsed optical states endowed with sizeable numbers of photons repre- 
sent a useful resource as they are robust with respect to losses and their reliable experimental 
detection and characterization is relevant, especially for establishing nonclassicality. 

In this Letter we report on a direct detection scheme aimed at measuring high-order 
correlations by means of a pair of hybrid photodetectors, which are photon-counting de- 
tectors endowed with a good linear response up to tens of photons [H]. In particular, we 
show that correlated bipartite optical states, namely a multimode twin-beam state (TWB) 
and a bipartite pseudothermal state generated in the mesoscopic photon-number domain 
[15] . can be effectively detected and characterized. To this aim, we define and derive the 
analytical expression of correlation functions at any order by only using quantities that can 
be experimentally accessed by direct detection, taking into account the non-unit quantum 
efficiency of the detection scheme. We also show that high-order correlations represent a 
useful discriminating tool of the nature of the state and demonstrate that, at increasing 
correlation order, the differences between classical and quantum states become more and 
more evident. 

The correlation functions are usually defined in terms of the normally-ordered creation 
and annihilation operators [16J: 
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where a& is the field operator of the mode fc-th and = ayi^. have a well-recognized 
meaning in connection with coherence properties of light and n-photon absorption process 
|17j . However, in a realistic direct detection scheme, we have only access to the shot-by- 
shot detected photons. We thus need to build analogous correlation functions for directly 
measurable quantities, namely: 

jk _ / X 

where rhk is the operator describing the actual number of detected photons in the fc-th 
arm of the bipartite state. Indeed, the functions g-£ can be linked to the g^ of Eq. ([!]), 
provided that a suitable description of the actual operation performed by the detector is 
given. If we assimilate the real detection to a Bernoullian process having efficiency rj, we 
can express all the operatorial moments of the detected-photon distribution as a function of 
those of the photon distribution, i.e., rh| = Ylh=i c ft(^)^fc' where the coefficients 0^(77) are 
given in Ref. [18]. Thus, by inserting this expression in Eq. ^ and exploiting commutation 
rules, we can express detected-photon correlations in terms of normally-ordered ones, i.e., 
9ih = Est=o £s J 9h-> where = 1 and e s>t depend on the physical parameters of the system 
under investigation. 

On the other hand, the expected results for correlations should be evaluated by taking 
into account all the realistic experimental conditions, that is not only imperfect detection 
and imbalance of the arms of the bipartite state, but also its multimode nature. Thus, we 
need to extend the theoretical description of the detection process to the case of multimode 
states, for which all the /i modes in the field are measured shot- by-shot. 

We consider a multimode TWB = <8)£ =1 \ip)k i n which each of the \i modes is in the 
same state, i.e., \ip)k = S n WV(l + (^)) n+1 |^)fc<8> \n)k, Vfc, and, thus, is equally populated 
with the same average number of photons per mode (n) = (n^), VA; [1J. By exploiting the 
pairwise correlations, we can state that the overall number of photons in each shot is the same 
in the two arms, being the sum of \i equal contributions coming from the /1 modes. We can 
thus write the multimode TWB in the following compact form [i/j^ } = Yl^=o VP" \ n ® ) ® \ n ® ) 
, where |n®) = S(n — Y2h=i n h) ®k=i \ n k)k represents the overall n photons coming from the 
H modes that impinge on the detector and: 

Pn n!Oi-l)!(^>/A*+irW<"> + l) B 
is the photon- number probability distribution for the multimode TWB. As one may expect, 
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we obtain a final result that depends only on the number of modes /x, the mean value of the 
number of photons (n) , and the overall detection efficiencies of the two detection chains, 771 
and r} 2 - In the case of 771 = r] 2 and considering correlation up to the 4-th order, we have (we 
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where (771) is the average number of detected photons and G k = Yl^iti + f 1 )/^- Indeed, the 
previous equations also describes two other classes of states: by choosing 77 = we obtain 
the correlation functions for a multimode thermal state, while for 77 = and fi — > 00 we get 
the results for the coherent state 
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FIG. 1: Scheme of the experimental setup. See the text for details. 

We generated a multimode TWB by the third harmonics (349 nm wavelength) of a mode- 
locked Nd:YLF laser regeneratively amplified at 500 Hz (High-Q Laser Production, Austria) 
impinging on a type-I BBO crystal (/3-BaB 2 04, Castech, China, cut angle 34 deg, 4-mm 
thick). According to the experimental setup sketched in Fig. [TJ we adopted a non-collinear 
interaction geometry in order to avoid possible residues of the pump beam. To obtain a good 
balancing of the quantum efficiencies, we selected two portions of the signal and idler cones 
close to frequency degeneracy, namely at 690 nm and 706 nm, by using two interference filters 
(IF, in Fig. [I]). The collection of a single coherence area in the two parties of the TWB state 
was obtained by inserting two pin-holes (PH, 2 mm diameter) at 107 cm and 109.5 cm from 
the BBO, respectively. The light passing the pin-holes was delivered through two multimode 
optical fibers to two hybrid photodetectors (HPD, R10467U-40, Hamamatsu, Japan), which 
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are detectors endowed with a partial photon-counting capability and a good linear response 
up to 100 photons. Their outputs were amplified (preamplifier A250 plus amplifier A275, 
Amptek), synchronously integrated (SGI, SR250, Stanford) and digitized (AT-MIO-16E-1, 
National Instruments). Each experimental run was performed on 50 000 subsequent laser 
shots at different values of the pump intensity. The self-consistent procedure to analyze the 
outputs of each detection chain is explained in details in Refs. [HI |2T] , here we only remind 
that the procedure allows us to obtain the mean value of detected photons (m), the number 
of modes /z, and the quantum efficiency rj directly from the experimental data. 

In previous works, we have already demonstrated that with the present detection appa- 
ratus assisted by the self-consistent analysis method it is possible to reconstruct both the 
detected-photon statistics [22] and the shot-by-shot second-order correlation in the number 
of detected photons [23] . Moreover, we have shown that the nonclassical nature of the TWB 
state can be proven by evaluating the noise reduction factor [15] . Nevertheless, as our TWB 
state is intrinsically multimode (fi ~ 100) and the effective detection efficiency is rather low 
(rj pa 5 %), the experimental characterization of the state in terms of sub-shot-noise correla- 
tions becomes challenging. On the contrary, the calculation of high-order correlations offers 
the possibility to better discriminate the nature of the state under investigation also in crit- 
ical situations. In Fig. [2] we plot the experimental data (full circles) obtained by evaluating 
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FIG. 2: High-order correlation functions for a multimode TWB state as a function of (mi + 1112) /2. 
Full circles: experimental data; Line+open symbols: theoretical expectations; Line+symbol: the- 
oretical prediction for a multimode thermal state. 

gi k up to the fourth order together with the theoretical expectations (line+open circles) for 
a multimode TWB calculated according to Eqs. Q by using the values of (m), [i and t] 
as directly obtained from the experimental data [151 18J . Since the values (mi) and (A2) 
are slightly different, in Fig. [2] and in all the other plots we use the quantity (mi + m 2 )/2. 
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The agreement of experimental data with theory is very high. For comparison, in Fig. [2] 
we also plot the theoretical predictions for a classically-correlated multimode thermal state 
having the same mean value and number of modes (line+symbols) J2D] • It is apparent that 
quantum correlations are always higher than classical ones, even if the absolute value of the 
difference is rather low due to the large number of modes and to the low value of quantum 
efficiency Moreover, as the order of correlations increases, the difference becomes larger. 
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FIG. 3: High-order correlation functions for coherent state (a) and single-mode pseudothermal 
state (6) as a function of (mi + 7712) /2 and multimode pseudothermal state (c) as a function of the 
number of modes. Full circles: experimental data; Line+open symbols: theoretical expectations. 
In panels (6) and (c) theoretical predictions for a TWB state with same parameters measured with 
different quantum efficiencies (?? = 5 %, dashed line, r\ = 30 %, dash-dotted line and r\ = 90 %, full 
line). Note the logarithmic scales in the plots. 

In order to explore the capability of our system in different regimes, we also measured 
a coherent state and single and multimode pseudothermal states generated by a rotating 
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ground-glass plate and divided at a beam splitter (23J, [21] . In Fig. [3] we plot the measured 
high-order correlations for some of these states. In panels (a) and (b), we plot the exper- 
imental results (full circles) and the theoretical expectations (line+open symbols) for the 
coherent state and for a single-mode thermal state as a function of the mean value of the 
number of detected photons on the two components of the bipartite state. In panel (c) 
we plot the experimental results (full circles) and the theoretical expectations (line+open 
circles) for a multimode thermal state at fixed mean value (rhi + 777,2) /2 = 3.88 as a function 
of the number of modes \i. All the experimental data superimpose to the theoretical expec- 
tations. As a comparison, in panels (b) and (c) we plot the value expected for a TWB state 
having the same physical parameters and measured with different quantum efficiencies. We 
note that the expected correlations for TWB are always higher than for classically corre- 
lated pseudothermal states, and the difference increases at increasing quantum efficiency. 
At fixed quantum efficiency, the maximum difference is achieved at any order for low values 
of both (m) and fi. We can conclude that by measuring high-order correlations we can 
achieve a better discrimination between quantum and classical correlations even at a very 
low quantum efficiency. 

A more direct way to establish nonclassicality is to address a suitable parameters satisfy- 
ing a boundary condition [25T427] . In our experiment, as a first check, we verify the Schwarz 
inequality: 

(mim 2 )/y '(mf) (ml) > 1, (5) 

the noise reduction factor: 

[{(mi - m 2 ) 2 ) - (rhi - m 2 } 2 ] / (mi + m 2 ) < 1 , (6) 
and the new inequality: 

22 _ r 131 r 121 

(mi)(m 2 ) 9 f Js + ^/HS> > 1 , (7) 

in which we have introduced symmetrized quantities [g ] s = \{g hk + g kh ) to take into 
account all the unavoidable asymmetries of the physical system. Indeed, all the above 
inequalities must be fulfilled by nonclassical light. Note that Ineq. ^ is equivalent to 
require g? 2 — g? 1 > in the case of perfectly balanced efficiencies. In Fig. [4] we plot the 
results for the three inequalities of above for our experimental data on multimode TWB. 
We note that, as expected, the measured state is nonclassical according to all the three 
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FIG. 4: Nonclassicality criteria for a multimode TWB as a function of the mean value of the number 
of detected photons. Full circles: experimental data; Line+open symbol: theoretical prediction. 

inequalities, but the amount of nonclassicality, that we can quantify with respect to the 
boundary unitary value, is much larger for the condition in Ineq. ([7]). This indicates that 
using high-order correlations can be helpful in all the situations in which the nonclassical 
nature is critical to be proven by other criteria, such as the case of a TWB state endowed 
with a high number of modes and/or a low quantum efficiency, for which we expect very 
small values of noise reduction factor. 

In conclusion, we have defined correlation functions at any order by means of quantities 
that can be experimentally accessed by direct detection with an experimental apparatus 
that allows accessing sizeable mean photon values. We have introduced a nonclassicality 
criterion based on these high-order correlations that represents a useful discriminating tool 
of the nature of the state in critical cases, in which other criteria are violated by a small 
amount: in some sense this criterion acts as an "amplifier" of nonclassicality violation. 
We have tested our theoretical and experimental procedure on a multimode TWB state by 
developing a multimode description that makes the calculation of high-order correlations 
straightforward. For the sake of completeness, we have also presented the direct comparison 
with an equally populated multimode pseudo-thermal state. The experimental results are 
in very good agreement with the theory, thus encouraging the exploitation of our setup for 
the reliable experimental characterization of quantum states for quantum technology. 
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